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ABSTRACT 

A set of nonhomogeneous radiation arid 
outflow boundary conditions which automat- 
ically generate prescribed incoming acoustic 
or vorticity waves and, at the same time, 
are transparent to outgoing sound waves pro- 
duced internally in a finite computation do- 
proposed. This type of boundary 
condition is needed for the numerical solu- 
tion of many exterior aeroacoustics problems. 
In computational aeroacoustics, the compu- 
tation scheme must be as nondispersive and 
nondissipati ve as possible. It must also sup- 
port waves with wave speeds which are nearly 
the same as those of the original linearized 
Euler equations. To meet these requirements, 
a high-order/large-stencil scheme is necessary. 
The proposed nonhomogeneous radiation and 
outflow boundary conditions are designed pri- 
marily for use in conjunction with such high- 
order/large-stencil finite difference schemes. 

1. INTRODUCTION 

Many exterior aeroacoustics problems in- 
volve incoming acoustic or vorticity waves in- 
teracting with an aircraft engine or the body 
of the aircraft. An example of this type of 
problem that is of current interest is the noise 
generation by the ingestion of free stream tur- 
bulence into a fan engine. Another example is 
the scattering and shielding of sound waves by 
raft wings and fuselage. To simulate these 
problems using computational aeroacoustics 
methods, the incoming acoustic and vortic- 
ity waves must be generated by the boundary 
conditions imposed at the outer boundaries of 
the finite computation domain. In this paper, 
a set of nonhomogeneous radiation and out- 
flow boundary conditions, which, when used 
in conjunction with a high-order finite differ- 
ence scheme, automatically generates the de- 
sired incoming acoustic and vorticity waves is 
proposed. 


In the presence of a uniform mean flow, 
the linearized Euler equations support three 
independent types of small amplitude distur- 
bances. They are the acoustic, the vorticity 
and the entropy waves. These waves, to lin- 
ear order, are uncoupled and propagate with 
different characteristics and wave speeds (see 
e g., Tam and Webb (1993)). In the compu- 
tational aeroacoustics literature, there seems 
to be an absence of suggestions as to how to 
generate such disturbances in the form of in- 
coming waves. There are two intrinsic difficul- 
ties. First, the imposed boundary conditions 
must not only generate the prescribed incom- 
ing waves, they must also be transparent to 
outgoing acoustic disturbances produced in- 
side the computation domain. Second, be- 
cause the Euler equations support both acous- 
tic and vorticity waves at the same spatial lo- 
cation, it is extremely difficult to generate a 
single type of disturbance without also pro- 
ducing the other type. Such coupling usually 
occurs at the corner regions of the computa- 
tion domain. The proposed nonhomogeneous 
radiation and outflow boundary conditions are 
designed specifically to overcome these diffi- 
culties. 

In this work, we will use the Dispersion- 
Relation- Preserving (DRP) scheme of Tam 
and Webb (1993) for numerical computation. 
There are two reasons for choosing the DRP 
scheme. First, the DRP scheme was designed 
so that the dispersion relations of the finite 
difference equations are always formally iden- 
tical to those of the original partial differential 
equations. This not only makes the waves sup- 
ported by the numerical scheme an excellent 
approximation of those of the partial differen- 
tial equations, it also assures that there will be 
no wave mode coupling in the numerical sim- 
ulation. The second reason for choosing the 
DRP scheme is that the exact plane acoustic, 
vorticity and entropy wave solutions of the fi- 
nite difference equations can be found analyt- 
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icallv. T!iis is discussed in .Section 2 of litis 
paper. These exact solutions are used in for- 
mulating the nonhoinogeneous radiation and 
outflow boundary conditions. 

Over the years, there has been a num- 
ber of investigations devoted to the devel- 
opment of radiation/nonreflecting boundary 
conditions but without incoming waves. Re- 
cently, Hixon. Shih and Maukbadi (1995) per- 
formed a detailed evaluation of the suitabil- 
ity of using the Thompson (1987, 1990) quasi- 
one-dimensionaJ characteristic boundary con- 
dition, the Giles (1990) Fourier mode decom- 
position boundary treatment and the asymp- 
totic boundary conditions of Tam and Webb 
(1993) for computational aeroacoustics appli- 
cations. They concluded that for their test 
problem, the radiation and outflow boundary 
conditions of Tam and Webb gave the least 
reflections and provided the only acceptable 
boundary treatment. In this work, we fol- 
low their recommendation and use the asymp- 
totic boundary conditions of Tam and Webb 
to develop the nonhomogeneous radiation and 
outflow boundary conditions that would au- 
tomatically produce very accurate incoming 
acoustic and vorticity waves in a finite com- 
putation domain. 


2. ACOUSTIC AND VORTICITY WAVES 
ON A GRID 

We will consider small amplitude distur- 
bances superimposed on a uniform mean flow 
of Mach number 'Af as shown in Figure L. 
Such disturbances are governed by the lin- 
earized Euler equations. We will use dimen- 
sionless variables with length scale s Ax = 
Ay (the mesh spacing), velocity scale = a 0 
(speed of sound), time scale = pressure 
scale s P 0125 (where po is the gas density). 
The dimensionless linearized Euler equations 
are 


df SF nr 

Jr +R e 7 +s dt = 0 


where F, the solution vector, and matrices R 
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fn (2.2) A/ r and A/ y are the mean flow Mach 
numbers in the x- and y-directions. If o is the 
direction of the mean flow measured from the 
x-axis. then A f t - A f cos A f v = M sin o. 

Upon discretizing ( 2 . 1 ) according to the 
DRP scheme, the finite difference equations 
may be written in the following form 

= C +•*!>' C J) (m» 

1=0 

Subscripts (f, m) axe the spatial indices of the 
mesh in the x- and p-directions. Superscript 
n is the time level. At is the time step. The 
vector is given by 

j=~3 3 

+ (2-4) 

j=-3 

The last sum on the right side of (2.4) is the 
artificial damping terms. R± x is the mesh 
Reynolds number. The numerical values of 
the stencil coefficients aj, b } and d } may be 
found in Tam (1995). 

We will look for plane wave solutions of 
finite difference equations ( 2 . 3 ) and ( 2 . 4 ) over 
an infinite mesh. The analytical form of the 
finite difference solution is 

= Rc {f e-'i •*♦*»—"*« J . (2.5) 

In (2.5). Re{ } is the real part of. F is the 
wave amplitude vector (a constant vector). 
(a,J) are the wave numbers, is the angular 
frequency. Substitution of (2.5) into ( 2 . 3 ) and 
(2.4). after some algebra, the following matrix 
equation for F is derived. 



The matrix A is equal to 


The corresponding eigenvector is 
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The quantities (a,/?) and -J are the wavenum- 
bers and angular frequency of the computa- 
tion scheme (the DRP scheme). 0(a) and 
D{$) are the artificial selective damping func- 
tions (see Tam (1995)). 

For nontrivial solution of equation (2.6), 
the determinant of the coefficient matrix A 
must be equal to zero. This condition leads to 
the general dispersion relation 


Suppose the direction of propagation of the 
incoming plane acoustic waves is \ (measured 
from the x-axis) as shown in Figure 1. Here 
the direction of propagation is taken to be nor- 
mal to the lines of constant phase. It follows, 
therefore, that cr and p must be such that 

tan \ (2.12) 

a 

A plane acoustic wave may be character- 
ized by its frequency and the direction of 
propagation The wavenumbers (a.J) of 
such a wave can be found by solving equations 
(2.10) and (2.12) simultaneously. With the in- 
clusion of artificial damping in the finite dif- 
ference scheme, equation (2.4), the wavenum- 
bers (cr, J) of the wave are complex. Since 
only a small amount of damping is added, the 
imaginary parts are small. The wavenumbers 
of the corresponding acoustic wave solution of 
the linearized Euler equations are, however, 
real. 

2.2 Vor deity waves 

The dispersion relation for the vorticity 
waves of the finite difference scheme is given 
by setting Z to zero in (2.9); i.e.. 


Z(u' — a 2 — p ) — 0. (2.9) 

We wish to potnt out that the Dispersion- 
Relation- Preserving scheme assures that dis- 
persion relation (2.9) is the same as that of 
the waves of the original linearized Euler equa- 
tions provided the following substitutions are 
made: u — w, cr — or, p — p and V 0. 


2.1 Acoustic waves 

For acoustic waves, the dispersion rela- 
tion is obtained by equating the second factor 
of (2.D) to zero. This gives. 


u -\f r a - A/, P + 


-a 3 - J' = 0. 


( 2 . 10 ) 
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+ D(*J)) = 0. 
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The corresponding eigenfunction is 
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Let r be the angle between the lines of con- 
stant phase of the vorticity waves and the x- 
axis (see Figure t) then 


tan r — — cot \ = — 


a 

P 


(2.15) 
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A vorticily wave may he specified by its 
frequency «j and the phase angle i". The 
wavenumbers of such waves are found by solv- 
ing equations (2.13) and (2.15) simultane- 
ously. Because of the inclusion of artificial 
damping terms, the wavenumbers are again 
complex. 


3. INHOMOGENEOUS RADIATION AND 
OUTFLOW BOUNOARY CONDITIONS 

We will now consider how to carry out 
numerical simulation of exterior aeroacoustics 
problems involving incoming acoustic or vor- 
ticily waves. By necessity, the computation 
domain is finite. Without loss of generality, 
we will consider the flow configuration to be 
as shown in Figure 1. With the inflow as 
prescribed, the bottom and the left side of 
the computation domain constitute the inflow 
boundaries. The right and top sides are the 
outflow boundaries. 

In the absence of incoming waves, the 
radiation boundary conditions of Tam and 
Web (1993) may be imposed along the inflow 
boundaries. In polar coordinates ( r,8 ) with 
origin at the center of the computation do- 
main. these radiation boundary conditions for 
the outgoing waves, F 0 ut. may be written as. 


= -V'(*) 


(*♦£)■ 


The acoustic wave propagation velocity V{8) 
of (3.1 ) is equal to, 


V{$) = \I cos(0-d) 


+ [1 - A/" sin 3 (0 — o)] ♦ . 


Along the outflow boundaries. Tam and Webb 
showed that the outgoing disturbances may be 
comprised of acoustic, vorticity and entropy 
waves. By means of the asymptotic solutions 
of these waves, they demonstrated that for 
pressure, p ou t. the outflow boundary condi- 
tion is the same as (3.1). For the velocity com- 
ponents. (Uout. c 0 ut). the approximate outflow 
boundary conditions are. 


^ u oul 
dt ^out 
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3.1 Incoming acoustic waves 

As discussed in Seel ion 2, an incoming 
plane acoustic wave on a grid has a mathe- 
matical solution in the form 


where F. the eigenvector, is given by (2.11) 
and tlie wavenumbers (n.d) are to be found 
by solving equations (2.10) and (2.12) simul- 


taneously. Let F‘j"^ be the numerical solution 


of the discretized linearized Euler equations. 
At the inflow boundaries, F^ is made up of 
the incoming and the outgoing acoustic waves. 
By subtracting F of {3.-4) from F^. the out- 
going acoustic wave solution. F out , is 


= C- F - 


Now F oul , being purely outgoing waves, 
must satisfy radiation boundary condition 

(3.1) . It is to be noted that F and F^J, 

and hence F oul are defined only on the so- 
lution mesh. To implement (3.1), the spatial 
and time derivatives must be discretized first. 

(3.1) may be rewritten in the form 


— F , "‘ D ( p wndl )\ ] 

dt l <■"* V /I (3.G) 


where W is 


W = -V{8) co$0-^— -f sin^-“ + ~ 
oz oy 2r 


-I- damping terms. 


The discretized form of (3.6) following the 
DRP scheme is. 




+ Re< F (^.g) 










lit the inflow boundary region (see Figure 
'2); i.e.. the bottom three rows ami the left- 
most three columns of mesh points, backward 
difference stencils are to be used to approxi- 
mate the spatial derivatives of (3.7). The par- 
ticular backward difference stencils to be used 
depend on the location of the point relative to 
the boundary of the computation domain. For 
example, for the corner point '.4* with f — L 
and m = M as shown in Figure 2, the dis- 
cretized form of (3.7) is. 

cos(0i„u) X a ' j A 
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where (r£, , 0£,.\/) are the polar coordinates 

of the point \V. c^ 51 and </, are the coef- 
ficients of the 5-point and 3-poiut damping 
stencils (see Tam ( 1995}). Equation (3.8) sup- 
plemented by formulae for in the form of 

(3.9) is the desired nonhomogeneous radiation 
boundary conditions. 

On the outflow boundaries, a similar 
treatment as above, starting with radiation 
boundary condition (3.3) instead of (3.1), 
leads to a set of inhomogeneous outflow 
boundary conditions. 


To test whether nonhomogeneous radi- 
ation boundary conditions (3.8) and (3. it) 
and the corresponding outflow boundaryeon- 
ditions can. indeed, generate an accurate 
plane acoustic wave propagating across the 
finite compulation domain without, at the 
same time, generating spurious vorticitv or en- 
tropy waves, a series of numerical simulations 
have been carried out. In the simulations, a 
100 x 100 mesh was used. In the interior re- 
gion. time marching scheme (2.3) and (2.4) 
were used to time-stepping the solution to a 
periodic state. Along the inflow boundaries, 
the variables of the solution were updated us- 
ing nonhomogeneous radiation boundary con- 
ditions (3.8) and (3.9). Along the outflow 
boundaries, nonhomogeneous outflow bound- 
in' conditions were used. Figure 3 shows the 
time history of convergence to the exact finite 
difference acoustic wave solution for the case 
M = 0.5. 6 = 30 deg, \ = 60 deg. u = 1.2 
and jf— = 0.025. The ordinate is the max- 
imum error; defined as the absolute value of 
the maximum difference (over the entire com- 
putation domain) between the numerical solu- 
tion and the exact solution of the finite differ- 
ence equations. The abscissa is the number of 
time steps. For this simulation, a zero acous- 
tic wave initial condition was used. As can be 
seen, over time, the computed solution con- 
verges to the exact incoming wave solutiou to 
machine accuracy. 

3.2 Incoming vorticity waves 

A set of nonhomogeneous radiation and 
outflow boundary conditions for incoming vor- 
ticity waves can be formulated in exactly the 
same way as discussed above. 

4. APPLICATIONS 

In this section we report the results of 
applying the nonhomogeneous radiation and 
outflow boundary conditions developed in the 
previous section to the numerical simulation 
of the scattering of plane acoustic waves by a 
cvlinder and the generation of sound by the 
interaction of a flat plate in a gust. 

4 1 Scattering of acoustic waves by a solid cylin- 
der 

The problem of the scattering of acoustic 
waves by a solid cylinder was simulated nu- 
merically in the time domain. The incoming 
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acoustic ivav^ wrre generated by the non ho- 
mogeneous radiation boundary conditions of 
Section 3. A 320 x 320 mesh was used in the 
computation. The cylinder with a diameter of 
32 was placed in the center of the compulation 
domain. Plane acoustic waves with a wave 
length of 8 mesh spacings entered the com- 
putation domain from the left boundary. In 
the numerical simulation, the linearized Euler 
equations were solved using the DRP scheme. 
A Cartesian boundary treatment using ghost 
values of pressure to enforce the solid surface 
boundary condition of zero normal fluid veloc- 
ity developed recently by Kurbatskii and Tam 
(1996) was applied around the surface of the 
cylinder. A zero acoustic wave initial condi- 
tion was used to start the computation. The 
solution was marched in time until a time pe- 
riodic state was reached. 

Figure 4 shows the contours of zero pres- 
sure obtained by the numerical simulation. 
The zero pressure contours of the exact so- 
lution are also plotted in this figure as dotted 
curves. The agreement between the numerical 
and the exact solution is good so that the two 
sets of curves are almost indistinguishable. In 
Figure 4, the shadow zone behind the cylinder 
shows up prominently. There are strong scat- 
tered waves propagating backward and to the 
two sides of the cylinder. They are responsible 
for the wiggles of the zero pressure contours of 
the figure. The good agreement between the 
numerical and exact solution provides strong 
evidence that the proposed nonhomogeneous 
radiation boundary conditions are, indeed, ac- 
curate and effective. 

4 2 Sound generation by the interaction of a Bat 
plate in a gust 

As another applicatiou of the nonhomo- 
geneous radiation and outflow boundary con- 
ditions, numerical simulations of sound gen- 
eration by the interaction of a flat plate in a 
gust (vorticity waves) were performed. The 
case of .Vf = 0.5 and a vertical gust with 

U = 0. > = 0.1 sin [J (£-«)] 

as shown in Figure 5 was considered. A 
200 x 200 mesh was used. The flat plate 
had a length of 30. It was located on the 
x-axis in the center of the computation do- 
main. The incoming vorticity waves were con- 
verted into the computation domam by the 


mean flow from the left boundary. The DRP 
time marching scheme was again used. The 
wall boundary condition on the plate was en- 
forced by the ghost point method. Nonhomo- 
geneous radiation boundary conditions were 
imposed on the left, top and bottom boundary 
regions. Nonhomogeneous outflow boundary 
conditions were imposed on the right bound- 
ary. The numerical solution was time-stepped 
from a zero vorticity wave initial condition un- 
til a time periodic state was established. 

Figure 6 shows tiie computed directivity 
of acoustic radiation in the top half-plane. No 
exact solution of the present problem is avail- 
able for comparison. Shown in this figure also 
are the computed directivities using tivice and 
three times the original spatial resolution. Jt 
is clear from the numerical results that there 
is numerical convergence. This gust-plate in- 
teraction problem is a generic turbomachinery 
noise problem. Direct numerical simulation of 
this problem as a scattering problem in the 
time domain has never been carried out be- 
fore. 

5. SUMMARY 

A set of nonhomogeneous radiation and 
outflow boundary conditions designed for use 
in conjunction with high-order/large-stencil 
finite difference schemes has been developed. 
The nonhomogeneous boundary conditions 
generate the desired incoming acoustic or vor- 
ticity waves and. at the same time, allow 
the scattered or internally generated distur- 
bances to leave the computation domam with 
almost no reflections. These nonhomoge- 
neous boundary conditions have been applied 
successfully to two benchmark aeroacoustics 
problems. In this work, only two-dimensional 
problems have been considered. The method- 
ology, however, can easily be extended to- 
three-dimensional problems. 

ACKNOWLEDGMENT 

This work was supported by NASA Lan- 
gley Research Center Grant NAG 1-1776. 

REFERENCES 

Giles. M.B., 1990, “Non-Reflecting Boun- 
dary Conditions for Euler Equation Calcula- 
tions". AIAA J.. 2S, 2050-2058. 

Hixon, R.. Shih, S.H. and Mankbadi. 
R.R. , 1995. “Evaluation of Boundary Con- 


337 






Figure 4 . Map of the zero pressure con- 
tours at the beginning of a cycle associated 
with the scattering of plane acoustic waves by 
a solid cylinder. Wavelength = Diame- 
ter of the cylinder = 32 Aj. 








